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Abstract 

Synchronization control of coupled continuous-time linear systems is 
studied. For identical systems that are stabilizable, a linear feedback law 
obtained via algebraic Riccati equation is shown to synchronize any fixed 
directed network of any number of coupled systems provided that the 
coupling is strong enough. The strength of coupling is determined by 
the smallest distance of a nonzero eigenvalue of the coupling matrix to 
the imaginary axis. A dual problem where detectable systems that are 
coupled via their outputs is also considered and solved. 

1 Introduction 

Synchronization is, at one hand, a desired behaviour in many dynamical systems 
related to numerous technological applications [9l [10l [8]; and, at the other, a 
frequently-encountered phenomenon in biology [2TJ [25l [TJ . On top of that, it is 
an important system theoretical topic on its own right. For instance, there is 
nothing keeping us from seeing a simple Luenberger observer [12] with decaying 
error dynamics as a two-agent system where the agents globally synchronize. 
Since people lack anything but good reasons to investigate synchronization, a 
wealth of literature has been formed, mostly in recent times [6l [22], [26] . 

An essential problem from a control theory point of view is to find con- 
ditions that imply synchronization of a number of coupled individual systems. 
This problem, which is usually studied under the name synchronization stability, 
has been attacked by many and from various angles. Two cases are of partic- 
ular interest: (i) where the dynamics of individual systems are primitive (such 
as that of an integrator) yet the coupling between them is considered time- 
varying; and (ii) where the individual systems are let be more sophisticated, 
however, coupled via a fixed interconnection. The studies concentrated on the 
first case have resulted in the emergence of the area now known as consensus in 
multi-agent systems [HI [Ml [71 [Til HH1 03 H3] where fairly weak conditions on 
the interconnection have been established under which the states of individual 
systems converge to a common point that is fixed in space. The second case, 



into which the problem studied in this note falls, has also accommodated im- 
portant theoretical developments especially by using tools from algebraic graph 
theory 27J. Using Lyapunov functions, it has been shown that spectrum of the 
coupling matrix plays a crucial role in determining the stability of synchroniza- 
tion [221 [H] notwithstanding it need not necessarily be explicitly known [I] . It 
has also been shown that passivity theory can be useful in studying stability 
provided that the interconnection is symmetric [3] [17l [20] . 

In this note we consider identical individual system dynamics that are linear 
time-invariant x — Ax + Bu. Under the weakest possible assumption that 
pair (A, B) is stabilizable, we search for a feedback law k(A, B) which would 
guarantee asymptotic synchronization for any fixed (directed) interconnection 
of arbitrary number of coupled systems provided that the coupling is strong 
enough. Following the tradition, we use the spectral information of the coupling 
matrix to determine the strength of connectedness. Namely, the farther the 
second eigenvalue (with largest real part) from the imaginary axis the more 
connected the network. We show that a linear k solving a linear quadratic 
regulation (LQR) problem performs the task. 

It is worth noting that most of the existing work on synchronization focuses 
on analysis rather than design. Also, the individual system dynamics are usually 
taken to be stable (i.e. the trajectories of the uncoupled system are required to 
be bounded.) In those respects the issue we deal with in this work is relatively 
different. Of particular relevance to this note are the works [23] and [21]. The 
former provides, for linear individual system dynamics, a linear feedback law 
that guarantees synchronization for all connected interconnections (regardless 
of the strength of coupling) under the extra assumption that matrix A is neu- 
trally stable. The latter establishes sufficient conditions for nonlinear individual 
system dynamics so that synchronization is achieved for strong enough coupling. 

In the remainder of the paper we first provide notation and definitions. Then, 
in Section [3] we formalize the problem and state our objectives. In Section [4] 
we show via our main theorem that optimal control theory yields us a feedback 
law which serves our purpose, i.e. synchronizes coupled systems for all network 
topologies with strong enough coupling. Finally, in Section [5] a dual problem 
is formulated and solved. 

2 Notation and definitions 

Let K>o denote set of nonnegative real numbers and | • | 2-norm. For A 6 C let 
Re(A) denote the real part of A. Identity matrix in R™ x ™ is denoted by /„ and 
zero matrix in R mx ™ by mXn . Conjugate transpose of a matrix A is denoted by 
A H . Matrix A £ C nxn is Hurwitz if all of its eigenvalues have strictly negative 
real partsH Given B G R" xm , C G M mx ", and A £ M" x "; pair (A, B) is 
stabilizable if there exists K G R mx ™ such that A — BK is Hurwitz; pair (C, A) 

1 Note that A is Hurwitz if there exists a symmetric positive definite matrix P such that 
A H P + PA < 0. 
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is detectable if (A T , C T ) is stabilizable. Let 1 6 MP denote the vector with all 
entries equal to one. 

Kronecker product of A G C mx ™ and B E C pxq is 



auB ■ ■ ■ a\ n B 



iB ■ ■ ■ a mn B 



Kronecker product comes with the property (A <S> B)(C ® D) = (AC) ® (BD) 
(provided that products AC and BD are allowed.) 

A (directed) graph is a pair (A/", A) where A/" is a nonempty finite set (of 
nodes) and A is a finite collection of pairs (arcs) (rii, nj) with rii, rij G A/". A 
pat/i from ri\ to is a sequence of nodes {ni, n2, . . . , n^} such that (n^, rii + i) 
is an arc for i G {1, 2, ...,£— 1}. A graph is connected if it has a node to which 
there exists a path from every other nodeH 

The graph of a matrix T := [7^] G M pxp is the pair (A/ - , .4) where Af = 
{ni, n2, . . . , n p } and (n^, n_y) G .4 iff 7^ > 0. Matrix T is said to be connected 
if it satisfies: 

(i) 7ij > for i ^ j; 

(ii) each row sum equals 0; 

(iii) its graph is connected. 

A connected T has an eigenvalue at A = with eigenvector 1, i.e. IT = 0, 
and all its other eigenvalues have real parts strictly negatived When we write 
Re(A2(r)) we mean the real part of a nonzero eigenvalue of T closest to the 
imaginary axis. 

Given maps & : IR> — > R n for i G {1, 2, . . . , p} and a map f : K> — » R", 
the elements of the set {£i(-) : i = 1, 2, . . . , p} are said to synchronize to £(•) 
if \£,i(t) ^ Z(t)\ —> as t — > 00 for all i. The elements of the set {£i(-) : i — 
1, 2, . . . , p\ are said to synchronize if they synchronize to some £(•). 



3 Problem statement 

3.1 Systems under study 

We consider p identical linear systems 

±i = Axi + Bm , i = l, 2, (1) 

where Xi G K™ is the state and Ui G M. m is the input of the ith system. Matrices 
A and are of proper dimensions. The solution of ith system at time t > is 

2 Note that this definition of connectedness for directed graphs is weaker than strong connec- 
tivity and stronger than weak connectivity. In 28 an equivalent condition (for connectedness) 
is given as that the graph contains a spanning tree. 

3 For the sake of completeness we provide a proof of this fact in Appendix. 
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denoted by Xi (t) . In this paper we consider the case where at each time instant 
(only) the following information 

p 

Zi = ^Hjixj-Xi) (2) 
3=1 

is available to ith system to determine an input value where 7^ are the entries 
of the matrix T G M. pxp describing the network topology. Nondiagonal entries of 
r are nonnegative and each row sums up to zero. That is, the coupling between 
systems is diffusive. 

3.2 Assumptions made 

We only make the following assumption on systems ([1]) which will henceforth 
hold. 

(Al) Pair (A, B) is stabilizable. 

3.3 Objectives 

We have two objectives in this paper. 

(01) Show that for each 5 > there exists a linear feedback law K G R mx ™ 
such that, for all p and connected T G W xp with — Re(A2(r)) > S, solutions 
of systems (HJ for Ui = Kzi, where Zi is as in ([2]), globally (i.e. for all initial 
conditions) synchronize. 

(02) Compute one such K. 

4 Main result 

For later use in this section we first borrow a well-known result from optimal 
control theory [T^: Given a stabilizable pair (A, B), where A g IR nx ™ and 
B G M. nxn \ the following algebraic Riccati equation 

A T P + PA + J„ - PBB T P = (3) 

has a (unique) solution P = P T > 0. One can rewrite ([3]) as 

(A - BB T P) T P + P(A - BB T P) + (/„ + PBB T P) = 

whence we infer that A - BB T P is Hurwitz. 

Lemma 1 Let A 6 R nx " and B 6 K" x ™ 1 satisfy ([3]) for some symmetric pos- 
itive definite P. Then for all a > 1 and ui G R matrix A — (a + joj)BB T P is 
Hurwitz. 
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Proof. Let e := a - 1 > 0. Write 

(A - (a + juj)BB T P) H P + P(A - (a + juj)BB T P) 

= (A — (a — juj)BB T P) T P + P(A - (a + juj)BB T P) 
= (A - (1 + e)BB T P) T P + P(A - (1 + e)BB T P) 
= (A- BB T P) T P + P(A - BB T P) - 2ePBB T P 
= -I n - (1 + 2e)PBB T P . 

Finally, observe that ([4]) is nothing but (complex) Lyapunov equation. 
Below is our main result. 



(4) 



Theorem 1 Consider systems ([T]). Let K := B T P where P is the solution to 
©. Given S > 0, for all p and connected T e R pxp with -Re(A 2 (r)) > S, 
solutions Xi(-) for i = 1, 2, . . . , p and m = max{l, S~ 1 }Kzi, where Zi is as in 
((2]), globally synchronize to 



x(t) := (r T ® e At ) 



x p (0) 



where r £ M. p satisfies r T T = and r T l = 1. 

Proof. Let Kg := max{l, S^ 1 }K. Combine |1| and @ to obtain 

p 



= Axi + BK S ^2 7y (a 

3=1 



Xi) . 



(5) 



Stack individual system states as x 
© as 



„T1T 



Then we can express 



x = (I p (g) A + T ® BK s )x . 

Now let Y £ C px ( p - 1 \ £ C( p -V* p , V £ 
C (p-l)x(p-l) be guch that 



(6) 

<p , and upper triangular A £ 



V = [1 Y] , V- 



r 

W 



and 



V- l TV 



A 



Note that the diagonal entries of A are nothing but the nonzero eigenvalues of 
r which we know have real parts no greater than —5. Engage the change of 
variables v := (V^ 1 ® i«)x and modify © first into 

v = (I p ® A + V^TV ® BK 5 )y 
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and then into 



A 

0(p— l)nxn 'p 1 







nx(p-l)n 

A + A ® BKj 



(7) 



Observe that l v -\ <8>A+ A®BKg is upper block triangular with (block) diagonal 
entries of the form A + XiBKg for i = 2, 3, . . . , p with Re(Ai) < — S. Lcmma[T] 
implies therefore that l v -\ ® A + A ® -B-fCs is Hurwitz. Thus ||7J) implies 



V(t) 

as t — > oo which yields 







nx (p— l)n 



^(p-l)nxn l)nx (p— l)n 



v(0) 







|x(i) - (lr 1 



-At 



x(0)| 



0. 



Hence the result. ■ 
By Theorem [1] we attain our objectives. In the next section we provide a 
dual result which may be more useful in certain applications. 



5 Dual problem 

Let p identical linear systems be 

Xi = A Xi + tii , yi = B Xi , i = 1, 2, 



(8) 



where £C{ € M n is the state, U{ S K™ is the input, and S M. m is the output of 
the ith system. Matrices A T and B T are of proper dimensions and they make 
a detectable pair (B T , A T ). Now consider the case where at each time instant 
the following information 



3 = 1 



(9) 



is available to ith system to determine an input value. Not surprisingly, the 
following result accrues. 

Theorem 2 Consider systems ©. Let L := PB where P is the solution to 
©. Given 5 > 0, for all p and connected T e W xp with -Re(A 2 (T)) > 6, 
solutions Xi(-) for i = 1, 2, . . . , p and Ui — max{l, 5~ 1 }Lzi, where Zi is as in 
(0), globally synchronize to 



x(t) := {r 1 ®e A l ) 
where r € R p satisfies r T T — and r T l = 1 



xi(0) 
asp(0) 



G 



6 Conclusion 



For identical (unstable) linear systems, we have shown that no more than sta- 
bilizability (detectability) is necessary for an LQR-based linear feedback law to 
exist under which coupled systems synchronize for all network topologies with 
strong enough coupling. We have considered directed networks and measured 
the strength of coupling via the real part of a nonzero eigenvalue (of the coupling 
matrix) closest to the imaginary axis. 

A Proof of a fact 

Fact 1 A connected T € M. pxp has an eigenvalue at A = with eigenvector 1 
and all its other eigenvalues have strictly negative real parts. 

Proof. First part of the statement directly follows from the definition. Since 
the entries of each row of T sum up to zero, it trivially follows that Tl = 0. 
Consider the dynamical system 

x = Tx 

where x G W with entries i, £l for i = 1, 2, . . . , p. Since Tl = wc can write 

v 

3=1 

By definition 7^ > for i ^ j. That brings us the following. For each r > 0, 
Xi(t) G [mxtii Xi(r), maxiXi(r)] for all i and t > t. In addition, since the graph 
of r is connected, interval [mhii Xi(t), maxiXi(t)] must shrink to a point as 
t — > 00. In other words, there exists x £ [min^ a;, (0), max^ 2^(0)] such that 

lim Xi{t) = x (10) 

t— >oc 

for all i (see |13|.) That readily implies by simple stability arguments that if 
A G C is an eigenvalue of T then Re(A) < 0. Also, for Re(A) = 0, the size of the 
associated Jordan block cannot be greater than unity. Possibility of a purely 
imaginary eigenvalue means sustaining oscillations and is therefore ruled out by 
(fT0|) . The only case left unconsidered is a second eigenvalue at the origin. That 
would imply, since it cannot be related to a Jordan block of size two or greater, 
that there exists a nonzero (eigen)vector deP such that Tv = and v ^ al 
for any a £ K. That contradicts pH]) . too. ■ 
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